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r*| , Abstract 

We study the elastic scattering of the lightest neutralino with a nucleus in the 
^ ■ framework of the minimal supersymmetric standard model (MSSM) with explicit 

flavor preserving CP violation, including the one-loop CP-violating neutral Higgs- 
boson mixing effects induced dominantly by the CP phases in the top and bottom 
(s)quark sectors. We construct the most general form of the effective Lagrangian for 
the neutralino-nucleus scattering in the limit of vanishing momentum transfers and 
then we perform a comprehensive analysis of the effects of the complex CP phases 
on the mass spectra of the lightest neutralino, neutral Higgs bosons and top squarks, 
and on the the spin-dependent and spin-independent neutralino-nucleus scattering 
cross section for three neucleus targets F, Si and Ge. The CP phases can reduce or 
enhance the neutralino-nucleus cross sections significantly, depending on the values 
of the real parameters in the MSSM. 
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1 Introduction 



A lot of cosmological and astrophysical observations have revealed that the dominant frac- 
tion of the matter in the universe and in the galactic halo is not luminous or dark [|]]. 
Furthermore, big-bang nucleosynthesis and galactic structure formation require some non- 
baryonic and non-relativistic matter for which various particle physics theories, in partic- 
ular, supersymmetric theories provide good candidates [El 0]. The lightest supersymmetric 
particle (LSP), which is guaranteed by R-parity conservation, is one of the most well- 
established and studied weakly interacting massive particle (WIMP) candidates of dark 
matter (DM). The null experimental results in searching for the anomalously heavy isotopes 
and the considerations in the most models suggest that the LSP in the minimal super- 
symmetric standard model (MSSM) might be the lightest neutralino, a linear combination 
of the four neutral superpartners of the U(l)y and SU(2) L gauge bosons and of the two 
neutral Higgs bosons; 

^ = N n B + N 12 W 3 + N 13 Hi + N u H° 2 , (1) 
where N is the unitary matrix diagonalizing the 4x4 neutralino mass matrix [|J. 

The lightest neutralino as a dominant cold dark matter (CDM) has very weak, but 
finite coupling strengths to quarks so that there might be a chance to detect the CDM 
neutralino even in laboratories directly or indirectly ||. If the couplings are too small, the 
neutralinos would not have annihilated in the early universe and it would be too much abun- 
dant today. Certainly, the capabilities of both direct and indirect search experiments for 
the lightest neutralino rely crucially on the size of the elastic neutralino-matter scattering 
cross section, in particular the neutralino-nucleus cross section because the detection rates 
for either direct or indirect searches are proportional to the cross sections. Since the first 
rough estimates were made by Goodman and Witten ||, the estimates of the neutralino- 
nucleus scattering cross section have been continuously and significantly improved [|7|, §]. 
Furthermore, there are many on-going and planned experiments with the goal of direct or 
indirect detections of the CDM neutralinos among which some experiments have already 
reported interesting but controversial results 0. Therefore, it is very important and timely 
to evaluate the cross section for neutralino-nucleus elastic scattering by including all the 
possibly dominant factors in determining the elastic scattering cross section as well as the 
neutralino-neutralino annihilation cross sections in the supersymmetric theories. 

In the present work we re-investigate the elastic scattering of the neutralino-nucleus 
scattering in the framework of the minimal supersymmetric standard model (MSSM) with 
R-parity focusing on all the dominant effects of flavor-preserving CP-violating complex 
phases on the cross section ||. The effects of the phase $ M of the higgsino mass parameter 
fi on the neutralino-nucleus scattering as well as the neutralino relic density have been 
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already studied by several works pD[ . However, there can exist additional effects through 
two different types of CP-violating phases^; one is the phase $1 of the U(l)y gaugino mass 
Mi and the other are the phases of the trilinear terms Af in the sfermion mass matrices. 
The latter phases can cause the CP-violating mixing among three neutral Higgs bosons at 
one-loop level so that the couplings of the Higgs bosons to the (s)particles are significantly 
modified [JTI, VI, 13 . We note that the phases of the trilinear parameters for the first and 
second generation sfermions are strongly constrained by the stringent experimental bounds 
on the electron and neutron electric dipole moments (EDMs) unless the sfermions are too 
heavy. On the other hand, since the third generation sfermion sectors give only two-loop 
suppressed effects on the electron and neutron EDMs when any generational mixing among 
sfermions is neglected, the CP violating phases involving the third generation sfermions 
can be very large |ll|]. In addition, the most dominant contributions to the CP-violating 
neutral Higgs boson mixing are from the third generation sfermion sectors because of their 
large Yukawa couplings flll, [T2|] . So, including the CP-violating neutral Higgs boson mixing 



due to the phases of the third generation trilinear parameters and the phase $1 as well as 
the phase <& M we will provide a comprehensive analysis of the dependence of the neutralino- 
nucleus elastic scattering cross section on those phases in the present work. 



The typical size of the momentum transfers exchanged in the neutralino-nucleus elas- 
tic scattering is of the order of 10 KeV, much smaller than the neutralino mass of the 
order of 100 GeV. For such tiny momentum transfers, the interactions of the lightest neu- 
tralino with spin-1/2 quarks can be described by two effective four-Fermi current x current 
terms; the spin-dependent axial-vector x axial- vector four-Fermi Lagrangian and the spin- 
independent scalar x scalar four-Fermi Lagrangian. We note that the spin-independent 
term becomes more significant for the scattering of the neutralino off nuclei with a large 
mass number because each spin-independent neutralino-nucleon scattering contributes co- 
herently to the overall spin-independent cross section. As a result, the spin-independent 
part dominates the spin-dependent part for large target nuclei. 

In the CP-noninvariant case, the scalar — pseudoscalar mixing among neutral Higgs 
bosons in the MSSM modifies the scalar couplings of each neutral Higgs bosons to fermions 
and neutralinos, leading to the significant changes of the spin-independent neutralino- 
nucleus cross section. It is, therefore, worthwhile to make a systematic investigation of the 
effects of all the relevant CP-violating phases on the neutralino-nucleus scattering cross 
section in the framework of the MSSM. For a systematic quantitative investigation, our 
analysis in the present work is based on a specific framework with the following assumptions: 

• The first and second generation sfermions are assumed to be very heavy so that they 
are decoupled from the theory. In this case, there are no constraints on the CP 

*Without any loss of generality, the SU(2)l gaugino mass M 2 can be taken to be real 
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violating phases from the neutron and EDMs. On the other hand, the annihilation 
of the neutralinos into tau pairs through the exchange of relatively light scalar tau 
leptons guarantees that the cosmological constraints on the DM densities be satisfied. 

• The explicit CP violation in the Higgs sector induced through the CP-violating ra- 
diative corrections from the top and bottom squark sectors is included. 

• Simultaneously, the effects of the induced CP-violating phase between Higgs doublets 
on the chargino and neutralino systems are explicitly included. 

• It is necessary to avoid the possible constraints from the so-called Barr-Zee-type dia- 
grams to the electron and neutron EDMs as well as from the Higgs search experiments 
at LEP. We take two values 3 and 10 for tan j3, the ratio of the vacuum expectation 
values of two neutral Higgs fields. 



The rest of the paper is organized as follows. In Sec. 0, we present a brief review of 
the CP-violating mixing among scalar and pseudo-scalar Higgs bosons in the MSSM Higgs 
sector on the basis of the recent work by Choi, Drees and Lee [T2|]. In addition, we discuss 



the neutralino mixing including the CP-violating phase induced from the MSSM Higgs 
sector. In Sec. [| we derive the complete analytic expressions for the spin-dependent and 
spin-independent effective four-Fermi Lagrangian relevant to our calculations explicitly, 
taking into account all the Higgs-boson, scalar-quark and Z-boson exchange diagrams in 
the limit of vanishing momentum transfers. Section |4] is devoted to a detailed numerical 
analysis for investigating the dependence of the neutralino-nucleus elastic scattering cross 
section for the three target nuclei F, Si, and Ge on the CP violating phases as well as 
other relevant real SUSY parameters such as tan/3 by taking a few typical parameter sets 
which effectively cover the whole parameter space. Finally, we summarize our findings and 
conclude in Sec. [5| 



2 CP— violating mixing in the MSSM 
2.1 Neutral Higgs— boson mixing 



In this section, we give a brief review of the calculation [[12] of the Higgs-boson mass matrix 
based on the full one-loop effective potential, valid for all values of the relevant third- 
generation soft-breaking parameters. 

The MSSM contains two Higgs doublets Hi, H 2 , with hypercharges Y{Hi) = —Y(H 2 ) = 
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—1/2. Here we are only interested in the neutral components, which we write as 

V2 
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where 0i.2 and are real scalar fields. The constant phase £ can be set to zero at tree 
level, but will in general become non-zero once loop corrections are included. The mass 
matrix of the neutral Higgs bosons can be computed from the effective potential fH}] 
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where we have allowed the soft breaking parameter m 
we have introduced the quantities 
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where the symbols g and g' stand for the SU(2)x and U(l)y gauge couplings, respectively. 
Q in Eq. (|3|) is the renormalization scale; the parameters of the tree-level potential, in 
particular the parameters m\, m\ and m\ 2 , are running masses, taken at scale Q. The 
potential ([|) is then independent of Q, up to two-loop corrections. 

The matrix M in Eq. (§) is the field-dependent mass matrix of all modes that couple 
to the Higgs bosons. The by far dominant contributions come from the third generation 
quarks and squarks. The (real) masses of the former are given by 
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where hb and ht are the bottom and top Yukawa couplings. The corresponding squark mass 
matrices can be written as 
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Here, H® and if^ are given by Eq. (0) while m 2 and m 2 are as in Eq. ([5]) and T> has been 
defined in Eq. (|]). In Eq. (H) m~, m~ and m~ are real soft breaking parameters, A b and 
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At are complex soft breaking parameters, and /x is the complex supersymmetric Higgs(ino) 
mass parameter. 



The mass matrix of the neutral Higgs bosons can now be computed from the matrix 
of second derivatives of the potential (|3]), where (after taking the derivatives) m 2 , m\ and 
m 2 2 sin(£ + #12) are determined by the stationarity conditions. The massless state G° = 
ai cos/3 — 02 sin (3 is the would-be Goldstone mode "eaten" by the longitudinal Z boson. We 
are thus left with a squared mass matrix M. 2 H for the three states a = a± sin (3 + a 2 cos (3, <pi 
and 02- This matrix is real and symmetric, i.e. it has 6 independent entries. The diagonal 
entry for a reads: 



M 2 
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\h t \ 2 m 2 
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and the CP-violating entries of the mass matrix, which mix a with <fii and 02 read: 
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where Aj and A^, which describe the amount of CP violation in the squark mass, read 
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and g(x,y) = 2 — [(a? + y)/(x — y)} \og(x/y). The definition of the mass squared ra\ and 
the dimensionless quantities X ttb , R t h and R' tb as well as the other CP-preserving entries of 
the mass matrix squared hA\ can be found in Ref. ||12|| . As noted earlier, the size of these 
CP-violating entries is controlled by A^ and A^. 



The real and symmetric 3x3 matrix M. 2 H can be diagonalized by an 3 x 3 rotation O; 
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with the increasing order of the three mass eigenvalues, m 2 Hl < m 2 H2 < m^, taken as 
a convention. Note that the loop-corrected neutral Higgs-boson sector is determined by 
fixing the values of various parameters; tua, /•*, At, A b , a renormalization scale Q, tan/3, 
and the soft-breaking third generation sfermion masses, Wq, m^, and m^. The radiatively 
induced phase £ is no more an independent parameter and it can be absorbed into the def- 
inition of the n parameter so that the physically meaningful CP phases in the Higgs sector 
are the phases of the re-phasing invariant combinations A t fie 1 ^ and Abfie^. This neutral 
Higgs-boson mixing changes the couplings of the Higgs fields to fermions, gauge bosons, 
and Higgs fields themselves so that the effects of CP violation in the Higgs sector can be 
probed through various processes (T5|. 



2.2 Top and bottom squark mixing 

On the other hand, the 2x2 Hermitian top and bottom squark mass matrices squared M. 2 
and M\ can be obtained after plugging the vacuum expectation values for the Higgs fields: 



Ml 
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— m b ^A b + /i* e~*t tan (3j m 2 ^ + m 2 — |m| cos 2/3 



These Hermitian mass matrices can be diagonalized by the unitary matrices, U$ and C/g, 

(tH(t)> 

respectively, where ti i2 and &i j2 are the mass eigenstates with their masses, 2 and , 
respectively. We note that the off-diagonal terms of each mass matrix squared are propor- 
tional to the corresponding fermion mass so that their contributions can be significant only 
for the third generation sfermions - top and bottom squarks. 



2.3 Neutralino mixing 

In general, the induced phase £ in Eq. fl2|) remains as a non-trivial physical phase and leads 
to a modification in the neutralino mass matrix; analytically, the induced phase plays a role 
of rotating the vacuum expectation value v 2 into u 2 e 1 ^ so that it modifies the neutralino mass 
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matrix describing the gauginos and higgsinos through electroweak symmetry breaking. In 
the weak-interaction basis (B, W 3 , H®, H®) the 4x4 symmetric, but complex neutralino 
mass matrix reads 
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| M^e** 1 


-m z cp s w 





M 2 
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-m z Cp s w 
m z op c w 




\ m z sps w e 
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m z sp s w e 
-m z sp c w e 





(13) 



where the phases of the U(l) gaugino mass Mi and the higgsino mass parameter /i, {<3>i, $ M } 
are explicitly given. The neutralino mass matrix is diagonalized in a symmetric way through 
a unitary matrix TV; 



x? 
xi 



N 



IB \ 

W 3 



(14) 



/ 



with the increasing order of the mass eigenvalues, m^o < m^o < m^o < Tn^o, as a convention. 
One can find that after an appropriate field redefinition there are two rephrasing-invariant 
phases $i and $ M + £ in the neutralino sector, where $i is the phase of the U(l) gaugino 
mass Mi and $ M the phase of the the higgsino mass parameter /i. 




Figure 1: The six mechanisms contributing to the neutralino-quark elastic scattering pro- 
cess Xi Q — > Xi Qi the spin-1 Z exchange, the three spin-0 neutral-Higgs-boson exchanges 
and the two squark q 12 exchanges. Here, the index k denotes 1,2 or 3. 
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3 Neutralino-nucleus elastic cross section 



3.1 Feynman rules 

In this section, we present all the Feynman rules in terms of mass eigenstates that are 
necessary for the neutralino-quark elastic scattering process (see Fig. 1). 



The interactions of the neutral gauge bosons Z to quarks are described by the same 
Lagrangian as given in the SM: 

Cz qq = — ^— q 7m \(Q q sw ~ T 3 ) Pl + Q q s 2 w Pr] q Z», (15) 

SwCw 

where Q q and T3 are the electric charge and the isospin of the quark, respectively. The 
interactions of the neutral Higgs bosons with quarks are described by the Lagrangian 

hd h 
C>H kqq = —j= d[0 2k + i Oi k S/375] dH k - -j=u[0?, k + i lk 0,375] u H k , (16) 

where the Yukawa couplings of down- and up-type quarks are given by 

7 em d em u 

hd = -7= , K = — = . (17) 

On the other hand, the SU(2) and U(l) gauginos, W 3 and B, have zero isospin and electric 
charges so that they do not couple to the Z boson. However, the Higgsinos H® and 
have non-trivial isospins so that the couplings of the Z boson to the lightest neutralinos, 
which are Majorana fermions, are described by the Lagrangian 

e 



\N 13 \ Z - N l4 



(xV75X?)^. (18) 



4s w c w 

where the Majorana nature of the neutralino is reflected in the fact that only the axial 
vector coupling but not the vector coupling is allowed. 



On the other hand, the interactions of the lightest neutralino with a pair of quark and 
squark come from both the gauge and Yukawa interactions. We consider the general flavor- 
diagonal squark mixing as well as the neutralino mixing and then we obtain the following 
Lagrangian 

g+(l-2)} + H.c, (19) 

where the coefficients B l ^ R (i — 1, 2) in the mass eigenstate basis expressed in terms of 
the squark mixing angle 6 q and phase (f> q read 

B\ L = cos 6 q A\ L + sin 6 q Af L , B 1 ^ = cos 6 q A L q R + e"^ sin q Af R , 

Bf = cos 9 q Af L - sin 9 q A L q L , B* R = cos 6 q A RR - e i<p " sin 8 q A LR , (20) 
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with the coefficients 
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for the up-type quarks it, c and t and 
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t\V #11 , 



(22) 



for the down-type quarks d, s and b. The chirality-preserving coefficients A^ L,RR originate 
from the gauge interactions, but the chirality-flipping coefficients A q R ' from the Yukawa 
interactions, respectively. 

Finally, the interactions of the neutral Higgs bosons to the lightest neutralinos involve 
both the neutralino mixing and the Higgs boson mixing. For the sake of notation we intro- 
duce the symbol G k defined in terms of the induced phase £, the neutralino diagonalization 
matrix N, and the Higgs diagonalization matrix O: 



G k = (JVia - t w Nu) U (Wis s p + N M cp e^) Q Xk + N 13 2k + N M 3k 



(23) 



With this abbreviation, the interaction Lagrangian for the coupling of the Higgs boson to 
the lightest neutralino pair is cast into a simple form: 



2s 



J2 Xi[Me(G k )+i%m(G k ) l5 }x°iH k 



W fc=l,2,3 



(24) 



where 3fte(Gfc) and Qm(Gk) denote the real and imaginary parts of the coefficient G k in 
Eq. (P3|), respectively. 

It is now straightforward to derive the coefficients A q and B q involving the effective four- 
fermion Lagrangian in the non-relativistic limit for the neutralino-quark elastic scattering: 



C 



off 



(25) 



The t-channel Z-exchange diagram contributes to the spin-dependent part, and the t— 
channel Higgs-exchange diagrams to the spin-independent part, while the s-channel squark- 
exchange diagrams contribute to both the spin-dependent and spin-independent parts. 
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After appropriate Fierz transformations, the coefficients A q and B q in the non-relativistic 
limit read: 
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(26) 



where Gp is the Fermi constant, and h q and T 3 are the Yukawa coupling and the third isospin 
component of the quark q, respectively. We have confirmed that in the CP-invariant the- 
ories, the expressions ( p6|) are consistent with those in J|. Note that the coefficient A q 
contains the Z-exchange contribution as well as the squark-exchange contributions, while 
the coefficient B q has the Higgs-exchange contributions as well as the squark-exchange 
contributions. The spin-dependent part is not suppressed by quark masses so that it can 
be large for the lightest neutralino of higgsino type, i.e, for large values \N\ 3 \ and |iVi4|. 
On the other hand, the spin-independent part from the B q terms is always proportional 
to quark masses and in particular the Higgs-exchange contributions are sizable only when 
the lightest neutralino is a well-balanced mixture of gaugino and higgsino states as seen 
from Eq. (^) and the size of the scalar-pseudoscalar mixing. Therefore, it might be naively 
expected that the small first and second generation quark masses will give rise to very small 
spin-independent cross sections. However, the spin-independent cross section can be sig- 
nificantly enhanced by coherent neutralino-quark scattering effects so that they dominate 
over the spin-dependent cross section for heavy nuclei. 



3.2 Neutralino— nucleus elastic scattering cross sections 

In order to obtain the scattering cross section of the lightest neutralino off heavy nuclei we 
need to know the detailed information on the configuration of the protons and neutrons 
inside each heavy nucleus and on that of quarks and gluons inside each proton and neutron. 
In the present work, we will not touch on this issue in detail, but simply take 19 F, 29 Si, 
and 73 Ge as three heavy nuclei and use the parametrizations of the form factors for each 
neutralino-nucleus scattering cross section as presented in Ref. 0. 

The total cross section for each elastic scattering process xl N — > Xi N with iV=F, Si, 
and Ge at zero momentum transfers can be divided into the axial-vector and scalar parts: 

a = cr A + a s . (27) 
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First of all, let us consider the spin-dependent part that is parametrized as follows: 

a A = —G 2 F J(J + l)m 2 r A 2 , (28) 

7T 

where m r is the neutralino-nucleus reduced mass m r = m^om^/(m^o + m^) and J the 
total angular momentum of the nucleus, the value of which is 1/2 for 19 F and 29 Si, 9/2 for 
73 Ge, respectively, and the spin- and target-dependent quantity A is given by 

A = j[a p (S p ) + a n (S n )], (29) 

where (S p ) = (N\S P \N) and (S p ) = (N\S n \N) and are the expectation values of the spin 
content of the proton and neutron group in the nucleus, respectively. 

These expectation values (S Ptn ) are different for each target nucleus and their explicit 
values for the target nuclei 19 F, 29 Si and 73 Ge are given in the shell model by 

(S P ) F =+0.415, (S n ) F =-0.047, 

(Sp)* = -0.002, (S n ) Si = +0.130, 

(S p ) Gc = +0.011, (S n ) Gc = +0.491. (30) 

On the other hand, the coefficients a p and a n are parametrized in terms of the quark spin 
contents of the proton and neutron, A p and A™, respectively, and the effective axial-vector 
couplings A q as 

q=u,d,s V^^F q=u,d,s V^^F 

In the present work, we take the factors A^ ,n (q = u, d, s) [16| to be : 

A p u = +0.77, A p d = -0.38, A^ = -0.09, 

A£ = -0.38, A n d = +0.77, AI? = -0.09. (32) 

More detailed information on the expectation (S p , n ) in Eq. (|30|) and the factors A^ and A" 
in Eq. (|3^) can be found in Ref. 0. 



Secondly, the spin-independent cross section to which each neutralino-proton or neutralino- 
neutron scattering contributes in a coherent manner is parameterized as 



,2 

l2 



4 Tfl 

a JH± [Z f p + ( A -Z)f n ]\ (33) 

7T 
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where Z and A denote the atomic number and mass number of the nucleus, respectively. 
In the limit of m x <C rriq and m q <C wig, the effective couplings of the lightest neutralino to 
protons and neutrons, f p and f n , are approximated to be 

— = e /r — + 4 /r e — , (34) 

Pi" q=u,d,s 1 q=c,b,t 1 

up to the lowest order in l/rrig, where the the parameters /^'J 1 are defined as 

r hy q=u,d,s 

In our present numerical analysis these quantities are taken to be 

fl =0.019, P =0.041, P =0.140, 

J-u L d l s 

r = 0.023, P = 0.034, r =0.U0. (36) 



as in Ref. [|I~7jj . 

We note in passing that some contributions from the other interactions such as so-called 
twist -2 operators are numerically small. 



4 Numerical results 

We are now ready to present some numerical results. It is known that loop-induced CP 
violation in the Higgs sector can only be large if both \fi\ and \A t \ ( or \Ab\ if tan/3 ^> 1) 
are sizable. In addition, for moderate values of tan j3 the contributions from the (s)bottom 
sector are still quite small so that our numerical results are not sensitive to and A^. In 
light of these points, we therefore choose 

\A t \ = \A b \ = 2rriQ = 2m & = 2m 5 = 1.0 TeV, (37) 

and take the phases of A t and A b to be equal. Numerically we find that for the large squark 
masses in Eq. ( ft7| ) the squark-exchange contributions to the neutraluno-nucleus scattering 
are very small compared to the other Z-boson and Higgs-exchange contributions. 

Since the main focus will be on the effects of the three CP-violating phases {$i, $ M , $a} 
on the neutralino-nucleus scattering cross section, we do not perform any parameter scan 
on all the real SUSY parameters in the present work. Nevertheless, it is crucial to cover the 



13 



overall parameter space effectively. So, except for the gaugino mass unification condition 
for the moduli of the SU(2) and U(l) gaugino masses [] 

| Mi | = - tan 2 9 W M 2 « 0.5 M 2 , (38) 
we consider three different sets for the gaugino and higgsino mass parameters: 



Gaugino 


M 2 


= 100 


GeV, 


v»\ 


= 1.0 


TeV, 


Higgsino 


M 2 


= 500 


GeV, 


H 


= 100 


GeV, 


Mixed 


M 2 


= \iA ~- 


= 150 


GeV. 







As will be shown later, the lightest neutralino is Bino-like in the gaugino case, while it is 
higgsino-like in the Higgsino case. 



Considering the Higgs search experiments and the constraints on the CP-phases of the 
top and bottom squark sectors from the electron and neutron EDMs at the two-loop level, 
we take two typical values of tan/3; tan/3 = 3 and 10. In addition, we take two values, 
150 and 500 GeV for the pseudoscalar mass parameter m^; In the former case, all the 
neutral Higgs bosons have masses less than 200 GeV, while in the latter case two neutral 
Higgs bosons are much heavier than the lightest Higgs boson. As a result, a significant 
CP-violating Higgs-boson mixing is expected only for = 150 GeV. Finally, we take the 
running masses m t {m t ) = 165 GeV and m&(mb) = 4.2 GeV as the top and bottom quark 
masses. We note in passing that the spin-independent cross section is rather strongly de- 
pendent on these quark masses WE 



The loop-induced phase £ between the vacuum expectation values of two neutral Higgs 
bosons can be adjusted to be zero by taking an appropriate renormalization scheme. With 
this adjustment and the universality assumption for the trilinear parameters, we can find 
that the Higgs-boson sector as well as the top and bottom squark sectors involves only the 
combination = $^ + $ M of the two rephasing-invariant phases $^ and $ M , while the 
light neutralino sector involves both $ x and $ At , but not 



4.1 Mass spectra 



Based on the above parameter sets as in Eq. (p9|) , let us first investigate the effects of 
the rephasing-invariant CP phases on the mass spectra for the the lightest neutralino, top 
squarks and Higgs bosons. Figure 2 shows the dependence of the lightest neutralino mass 
on the phase ^ for five different values {0°, 45°, 90°, 135°, 180°} for the phase $ x ; the 

^The phase $i of Mi is zero up to the one-loop level for a positive Mi in the minimal supergravity 
model. 
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upper left (right) frame is for tan/3 = 3 (10) in the gaugino case, the middle left (right) 
frame is for tan (3 — 3 (10) in the higgsino case, and the lower left (right) frame for tan (3 = 3 
(10) in the mixed case. The approximate form of the lightest neutralino mass is given by 

f Gaugino : \M X \- sin 2(3 cos ($i + $ M ), 

Higgsino : \fj,\ — ^f- (1 + sin 2(3) 

in the gaugino and higgsino cases, respectively. 
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Figure 2: The dependence of the lightest neutralino mass m^o on the phase $ M for five 
different values {0°, 45°, 90°, 135°, 180°} for the phase $i; the upper left (right) frame is 
for tan (3 = 3 (10) in the gaugino case, the middle left (right) frame is for tan (3 = 3 (10) 
in the higgsino case, and the lower left (right) frame for tan (3 = 3 (10) in the mixed case. 



Combined with the approximate expressions fl40|) , Figure § shows several interesting 
features about the sensitivity of the the neutralino mass to the phases $a and $^1 
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• The mass is mainly determined by the modulus of Mi and /i in the gaugino and 
higgsino cases, respectively, while in the mixed case the mass is affected by both Mi 
and \i. 

• As the phase-dependent parts are accompanied by a factor sin 2/? as shown in the 
expressions, the sensitivity of the neutralino mass to the phases is smaller for larger 
tan f3 in every case. 

• In the mixed case, the sensitivity of the neutralino mass to the phase $1 is larger 
when $ M = or 2 7r. 

Consequently, it is expected that only in the mixed and higgsino cases with relatively small 
| Mi | and \n\ the neutralino-nucleus scattering cross sections can be affected by the phases 
$ M and $i through the lightest neutralino mass rriyO. 
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Figure 3: The dependence of the light and heavy top squark masses on the phase 
for (a) \fi\ = 1 TeV and tan/? = 3, (b) |//| = 100 GeV and tan/? = 3, (c) \fi\ = 1 TeV and 
tan/3 = 10, (d) \fi\ = 100 GeV and tan/3 = 10 with the parameter set ft3l[ ) for the squark 
mass parameters and trilinear parameters. 
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For moderate values of tan (3 the loop-induced effects to the neutral Higgs boson sector 
is dominated by the top (s)squark contributions due to the largest top Yukawa coupling. 
In this light, we present in Fig. |3| the dependence of the top squark masses on the 
phase ^Afj, for (a) \fi\ = 1 TeV and tan/3 = 3, (b) \fi\ = 100 GeV and tan/3 = 3, (c) = 1 
TeV and tan/? = 10, (d) = 100 GeV and tan/3 = 10 with the parameter set (|3~7"D for 
the squark mass parameters and trilinear parameters. Analytically, we find that the phase 
dependence of the top squark masses is determined by the quantity: 



A 12 = m t 5sm(A t fi e i? ) cot (3 



(41) 



Figure |3] and Eq. fl41~|) clearly show that the top squark masses are strongly dependent on 
the phase only in the gaugino case with large \fj,\ and small tan/3. In the other cases, 
the top squark masses are almost insensitive to the phase <&Ai±- 
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Figure 4: The Higgs-boson mass spectrum with respect to the phase ^a^ in (a) the gaugino 
case and (b) the higgsino case with tan/3 = 3, and in (c) the gaugino case and (d) the 
higgsino case with tan/3 = 10; the pseudoscalar mass parameter vtia is set to be 150 GeV. 
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The mass spectrum of the neutral Higgs bosons is sensitive to the combination = 
&a + 3>/i- So, we exhibit in Fig. 3 the Higgs-boson mass spectrum with respect to the 
phase with uia = 150 GeV. The mass difference between the the lightest Higgs boson 
and the heavier Higgs boson is more significant for smaller tan (3. The reason is that 
for the pseudoscalar mass comparable to mz, the off-diagonal entries of the Higgs-boson 
mass matrix are essentially proportional to sin 2/3. On the other hand, the sensitivity of 
the Higgs-boson masses to the phase Qa^ is strongly suppressed for large ttia due to the 
suppressed mixing. Because of this feature we have shown the Higgs mass spectrum only 
for ttia = 150 GeV in Fig. 3. For moderate values of tan/3 the size of the one-loop induced 
CP-violating neutral Higgs-boson mixing is dictated by the factor 

A t - = — ; - 2 ; . 42 

t-2 tl 

Therefore, only when both \A\ and are large, the Higgs boson masses become sensitive 
to the phase This feature is clearly reflected in Fig. 3 for the Higgs mass spectrum; 

in the gaugino case with a large value of = 1 TeV, the Higgs boson masses depend 
strongly on the phase, while in the higgsino case with a small value of = 100 GeV, the 
Higgs boson masses are almost insensitive to the phase. Consequently, the large effects of 
the phase $An on the Higgs masses and mixing are expected for large \fi\ as well as small 
rriA and tan j3. 



4.2 Neutralino— nucleus scattering cross sections 

With the previous comprehensive investigations of the mass spectra for the lightest neu- 
tralino, top squarks and neutral Higgs bosons, let us study in this section the dependence 
of the neutralino-nucleus cross sections on the CP-violating phases {$1, $ M , &a} as well as 
the other real SUSY parameters. 

Firstly, we evaluate the elastic cross sections in the CP-invariant theories. We take 
the three phases to be vanishing; $a = $ M = $1 = and consider F, Si and Ge as three 
typical target nuclei. Table 1 shows the neutralino-nucleus elastic scattering cross sections; 
spin-independent and spin-dependent parts, as and a a, for the three different scenarios 
and two values of m^; the upper table is for tan/5 = 3 and the lower part for tan/3 = 10. 
Note that the size of the cross sections depends significantly on the scenarios and on the 
target nuclei as well as on the pseudoscalar mass parameter m^: 

• The cross sections are strongly suppressed in the gaugino case, but they are very 
much enhanced in the mixed and higgsino cases. This suppression feature is more 
predominant in the spin-dependent cross section, for the coupling of the Z boson to 
the neutralinos is strongly suppressed in the gaugino case, i.e. the absolute magnitudes 
of 7V 13 and /V 14 are extremely small. 
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• The spin-independent cross sections are much more enhanced for heavier nucleus and 
become dominant over the spin-dependent cross section as can be expected from the 
coherent scattering effects. 

• The spin-independent cross sections are strongly suppressed for large pseudoscalar 
mass ttia and large tan (3, especially in the gaugino case. 

The elastic scattering process Xi Q ~ * Xi Q i s related to the neutralino-neutralino annihila- 
tion process, XiXi ~^ QQi contributing to the relic density of the SUSY dark matter. We 
note that in most cases the gaugino scenario is favored by the estimates of the relic density. 



Table 1: The neutralino-nucleus elastic scattering cross sections; spin-independent and 
spin-dependent parts for the three different scenarios and two values of itla in the CP- 
invariant case, $ M = $1 = $a = 0. F, Si and Ge are considered as three different target 
nuclei; the upper table is for tan/3 = 3 and the lower part for tan/5 = 10. 



tan (5 — 3 

m A [GeV] 


F 

<7 8 [fb] <7 A [fb] 


Si 

<7 8 [fb] (7 A [fb] 


Ge 

<7 8 [fb] <7 A [fb] 


Gaugino 150 
500 


0.17 0.02 
0.15 0.02 


0.63 0.002 
0.57 0.002 


11.6 0.028 
10.6 0.028 


Higgsino 150 
500 


10.7 77.0 
11.3 77.0 


43.3 8.79 
45.8 8.79 


992 161 
1049 161 


Mixed 150 
500 


25.0 120 
26.2 120 


94.7 12.8 
99.1 12.8 


1769 191 
1852 191 



tan (3 


= 10 




F 




Si 


Ge 




m A [GeV] 


<7 S [fb] 


<r A [fb] 


(7s [fb] 


<7 A [fb] 


<7 S [fb] <7 A [fb] 


Gaugino 


150 


0.38 


0.02 


1.45 


0.002 


26.9 0.035 




500 


0.003 


0.02 


0.01 


0.002 


0.20 0.035 


Higgsino 


150 


59.7 


114 


244 


13.1 


5741 247 




500 


60.3 


114 


246 


13.1 


5800 247 


Mixed 


150 


160 


157 


622 


17.3 


12689 282 




500 


161 


157 


626 


17.3 


12783 282 
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(a) Gaugino, tanpS=3 



(c) Gaugino, tan(3=10 
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Figure 5: The ratio 7£ Gc [$ M ; $,4, $i] unt/i respect to the phase <& M /or the parameter set 



(31) and ttla = 150 GeV in (a) the gaugino case and (b) the higgsino case with tan/3 = 3 ; 
and in (c) the gaugino case and (d) the higgsino case with tan/3 = 10. Here, the phase $i 
is taken to be 0°. 



Secondly, we study the effects of the CP-violating phases on the neutralino-nucleus 
scattering cross sections. Note that the spin-dependent cross sections a A are independent of 
&a and almost independent of the phases <3>i and $ M , especially in the gaugino and higgsino 
cases, because the neutralino mixing elements iV 13 and iV 14 determining the coupling of the 
Z boson to the lightest neutralinos are suppressed by the factors mz/\n\ and m z /\Mi^\ in 
the gaugino and higgsino cases, respectively. Even in the mixed case the lightest neutralino 
is almost Bino-like because of the gaugino mass unification condition so that this scenario 
also causes the spin-dependent cross sections are almost independent of the CP-violating 
phases. Taking into account these points, we consider the ratio of the spin-independent 
cross sections only for the heaviest target nucleus Ge : 



o- s [^;$a,$i] 
a s [0 ; 0,0] 



(43) 
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for several values of the phases {$1, $,4} and present the numerical results for the spin- 
independent cross section ratio (|43| ) by taking rriA = 150 GeV and the parameter set (|37|) 
in three figures; the first figure is for $x = 0° (Figure 5), the second one for $! = 90° 
(Figure 6) and the third one for $ a = 180° (Figure 7), respectively. Each figure contains 
four frames for (a) the gaugino case and tan/3 = 3, (b) the higgsino case and tan/5 = 3, 
(c) the gaugino case and tan/3 = 10 and (d) the higgsino case and tan/3 = 10; in each 
frame the solid line is for &a = 0°, the dotted line for $,4 = 90° and the dashed line for 
$4 = 180°. 
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Figure 6: The same ratio 7?. Gc [$^5 $^,90°] as in Fig. I but for $ x = 90 c 



Comparing the results presented in the three figures, we find several interesting features 
about the dependence of the spin-independent ratio 7l Gc : 

• In all the figures and all the scenarios the spin-independent ratio is strongly dependent 
on the phase it is striking that certain values of $ M render the cross sections almost 
vanishing. 
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• The ratio is very sensitive to the phase §a in the gaugino case, where the Higgs-boson 
couplings to the lightest neutralinos as well as the quarks are significantly affected 
by the phases &a and $ M as well as the phase $1. This $^ dependence is expected 
to be more significant for small rriA due to a larger mixing between the scalar and 
pseudoscalar Higgs bosons. The ratio itself varies more significantly with the phases 
$ M and $a for tan /? = 10 than tan (3 = 3. 

• In the higgsino case, the cross section ratio IZce is (almost) always suppressed for 
non-trivial values of $ M , i.e. the neutralino-nucleus cross section is maximal near 
$ M = and 360°. 

• On the contrary, the cross section ratio is enhanced or suppressed depending on the 
values of the phases {$i, $ M , $,4} in the gaugino case. Interestingly, the cross section 
ratio is maxmal when the sum of the phases $ M and $1 is 0° (360°) or 180°, while 
the ratio is almost vanishing when the sum of the phases $ M and $1 is 90° or 270°. 
This feature can be understood from the the fact that in the gaugino case the matrix 
elements Nik = 1 - 4) as well as the lightest neutralino mass depend on only the 
rephasing-invariant combination ^ + $! with a good approximation. 

• Comparing the lines for three different values of §a in each frame of Figs. ||, || and 
in the gaugino case, we can confirm that the phase determining the CP-violating 
Higgs-boson mixing is indeed the rephasing-invariant combination §Afj, of the phases 
$ A and $ M . 

All these features can be understood by noting that the spin-independent interactions 
are primarily determined by the coupling strengths of the Higgs bosons to the neutralinos 
and the quarks given in Eqs. (|i~6| ) and (j2"4|). The scalar couplings of a Higgs bosons to quarks 
can vanish when the Higgs boson becomes (almost) a CP-odd state for certain values of 
&A.fj, through the CP-violating scalar-pseudoscalar mixing. Certainly, for some values of 
the phases the scalar couplings can be enhanced compared to the CP-invariant case. On 
the other hand, the scalar couplings of the Higgs bosons to the lightest neutralinos, which 
is given by the real part of Gk in Eq. (|23|) , involves not only the CP-violating Higgs- 
boson mixing but also the neutralino mixing. Therefore, these couplings are expected to 
be strongly dependent on all the phases in the gaugino case. 

To summarize, the CP-violating phases could reduce or enhance the neutralino-nucleus 
cross sections compared to those in the CP-invariant theories. The phase $ M as well as the 
phase $1 affects the cross section significantly in most cases, while the phase &a modifies 
the (spin-independent) cross section in the gaugino case with large through the phase 
combination = $^4 + $ M for relatively small itla- 
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Figure 7: The same ratio 7^ Gc [$ M ; $i] as m i*%s. [| and frirf /or $i = 180 c 



5 Summary and conclusions 

In this paper, we have made a comprehensive investigation of the effects of the CP-violating 
phases on the neutralino-nucleus elastic scattering cross sections in the framework of the 
MSSM through the neutralino mixing and the one-loop induced Higgs-boson mixing from 
mainly the top (s) quark sector due to the largest top Yukawa coupling. For the sake of 
numerical analysis, we have imposed the universal relation for the squark mass parameters 
and the gaugino mass unification condition only for the moduli of the SU(2) and U(l) gaug- 
ino masses. We have taken large sfermion masses and moderate values of tan/3, allowing 
the CP-violating phases to have large non-trivial values without violating the stringent 
constraints of the electron and neutron EDM measurements, as well as large trilinear pa- 
rameters for a significant neutral Higgs-boson mixing. 
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Based on three scenarios - gaugino, higgsino and mixed - for the neutralino sector, we 
have performed a detailed numerical analysis for the mass spectra of the lightest neutralino, 
top squarks and neutral Higgs bosons. The lightest neutralino mass is strongly dependent 
on the phases $1 and $ M for small tan/5 and in the mixed case with comparable gaugino 
and higgsino mass parameters. On the contrary, the top squark and Higgs boson masses 
are significantly affected by the combination of the phases $^ and &a in the gaugino 
case with large and with small tan/?. 

After estimating the neutralino-nucleus scattering cross sections for the target nuclei, 
F, Si and Ge for the three scenarios in Eq. (139) in the CP-invariant theories, we have 
studied the effects of the CP-violating phases on the spin-independent cross section ratio 
TZ Gc for the target nucleus Ge, for which the spin-independent cross section is significantly 
enhanced by large coherent scattering effects due to its large mass number. We have found 
that the CP-violating phases could reduce or enhance the cross section ratio significantly. 
The phase $ M as well as the phase $1 affects the ratio section in almost all the cases, while 
the phase &a modifies the (spin-independent) cross section in the gaugino case with large 
\n\ through the phase combination <&An = $^ + $ M , when the pseudoscalar mass tua is 
relatively small. 

Consequently, the CP-violating phases as well as all the real parameters in the MSSM 
are very crucial in estimating the neutralino-nucleus cross sections and so in determining 
the possibility of experimentally detecting the lightest neutralinos if they indeed constitute 
a major component of cold dark matter in the Universe. 
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